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In this paper, we propose a unified formalism, using Green's functions, to integrate out the 
electrons in an insulator under uniform electromagnetic fields. We derive a perturbative formula for 
the Green's function in the presence of uniform magnetic or electric fields. Applying the formula, 
we derive the formula for the polarization, the orbital magnetization, and the orbital magneto- 
polarizability, without assuming time reversal symmetry. Specifically, we realize that the terms linear 
in the electric field can only be expressed in terms of the Green's functions in one extra dimension. 
This observation directly leads to the result that the coefficient of the 6 term in any dimensions is 
given by a Wess-Zumino-Witten-type term, integrated in the extended space, interpolating between 
the original physical Brillouin zone and a trivial system, with the group element replaced by the 
Green's function. This generalizes an earlier result for the case of time reversal invariance [see Z. 
Wang, X.-L. Qi, and S.-C. Zhang, Phys. Rev. Lett. 105, 256803 (2010)]. 



I. INTRODUCTION 

In recent years, there has been a series of development 
on how to calculate the polarization and orbital mag- 
netization quantum-mechanically in insulating systems, 
in terms of Bloch wave functions<^"— The polarization P 
measures the position differences between the band elec- 
trons and the lattice ions, which is in general nonzero in 
a crystal without inversion symmetry. In an open sys- 
tem, it results in boundary charges, which gives rise to 
the energy density AE' = — P ■ E in an external elec- 
tric field. Inside the bulk, it is only defined modulo ea 
where e is the electron charge, and a is the lattice spac- 
ing. This arises naturally from the ambiguity of associ- 
ating an electron with a given ion. The reason why it has 
not been fully computed until rather recently is due to 
the unboundedness of the operator r, which is the nat- 
ural quantity to calculate the expectation value, when 
the lattice ions stay fixed. In Refi^, they overcome this 
problem utilizing the Wannier orbitals. The Wannier or- 
bitals are localized in an insulating system, and the ma- 
trix elements between them are well-defined, even if the 
operator itself is unbounded. They also provide an alter- 
native method, which is to compute the current response 
under a change of the Hamiltonian without breaking the 
lattice translation symmetry. The integration of the cur- 
rent gives the change of the polarization. This method 
avoids calculating (r) altogether, at the expense of intro- 
ducing an extra parameter and a Hamiltonian depending 
on it. This method also leaves the impression that only 
the diffrencc of the polarization is properly defined. We, 
however, have argued otherwise in RefJ^. 

Similarly, the magnetization M is generally nonzero 
when the time-reversal symmetry is absent. The energy 
density change in the presence of an external magnetic 
field is given by IS.E = — M • B. Usually the time-reversal 
symmetry is broken due to magnetism, and the domi- 
nant contribution to the magnetization is of spin origin. 
Nevertheless, via spin-orbit coupling the orbital motion 



contributes to the magnetization as well. The orbital 
magnetization is first computed in Refi^i^, either by semi- 
classical methods^, or by calculating the matrix elements 
of r X V using again the Wannier orbitals^ In the latter 
work, however, special care has to be taken toward the 
boundary, as they show that an extra term Mic arising 
from the "itinerent current" fiowing around the bound- 
ary has to be included, in addition to the "local current" 
contribution AfLCi which involves the matrix element of 
r X V between the bulk Wannier orbitals, to give the cor- 
rect and the same answer as in the former calculation. 
Later on in Ref.— , the authors give a full quantum me- 
chanical derivation, which calculates the energy of the 
system in a finite magnetic field based on the finite q 
perturbation theory of the vector potential, and taking 
q in the end. By taking a derivative with respect to 
B one gets the magnetization M. 

After the discovery of the three-dimensional (3D) topo- 
logical insulator—, it was shown in Refi^ via a dimen- 
sional reduction procedure from one higher dimension, 
that the 3D topological insulator is characterized by the 
6'-term, Cg = If^-E • B, with 6* = tt in the effective the- 
ory, where the electrons arc integrated out. With time 
reversal symmetry, the coefficient 9 is given by the inte- 
gral of the Chcrn-Simons term of the Berry's connection 
in the Brillouin zone, which is independent of the gap 
size, has a 27r ambiguity, and therefore can only take the 
value or tt. In Reffi, it is shown that with time rever- 
sal symmetry, the 6 coefficient can also be written as a 
Wess-Zumino-Witten- (WZW-) type term, integrated in 
an extended space, where one interpolates between the 
system in question and a trivial insulator. 

Since then, there have been various attempts to calcu- 
late this term explicitly in three spatial dimension o^°'^^ , 
without assuming time reversal symmetry. The 6-icvm. 
is treated as an energy density, which can be under- 
stood as either the polarization in a magnetic field or the 
magnetization in an electric field, which are denoted as 
orbital magneto-polarizability (OMP) or orbital electro- 
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susceptibility (OES). In Rcfi^i, they offer two ways to cal- 
culate the OMP, similar to the methods described above 
for the polarization. In the first method they again cal- 
culate the current response to a change of Hamiltonian, 
but with a small uniform magnetic field turned on all the 
time. To perturb in the magnetic field they use the den- 
sity matrix perturbation theory. In the second method, 
they evaluate the matrix element r between Wannier or- 
bitals in the magnetic field, using the finite q perturba- 
tion theory and then take the q — ^ limit. In general, 
however, they find that = a'-B^; that is, the OMP 
is not diagonal, and the O-term is just a part of it. This 
result is confirmed by RefiiS, where they use the Wannier 
orbitals to study the OES, by calculating the magnetiza- 
tion in an electric field. In this calculation, similar care 
has to be taken on the boundary. In either calculation, 
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the diagonal part of a* is defined modulo |^ . The easiest 
way to understand this physically is the following. If we 
consider a cylinder geometry with the material in ques- 
tion in the bulk, an integer quantum Hall layer on the 
surface of the cylinder will change the diagonal response 
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by 2||r, where n is the filling factor of the layer. 

While all the results in the end agree with each other, 
the derivations are diverse, with various limitations and 
subtleties, as noted below: 

(i) The Wannier orbitals can only be defined when the 
Chern number of the bands is zero.— Furthermore, to 
use the operator r or r x v, one is essentially limited to 
settings with open boundary conditions, as they are not 
well-defined on a torus. The boundary then has to be 
treated carefully, even if we are only interested in bulk 
properties: in Ref.— , to get the correct expression for the 
magnetization, they have to consider two contributions, 
Mlc and Mjc as we briefiy mentioned. The first term 
is the usual local matrix element between the bulk Wan- 
nier orbitals and the second term is itinerent, comes from 
the boundary, where the Wannier orbitals are deformed. 
The second term, however, can be written as a function 
of bulk parameters and is then argued to be present even 
in a setting with periodic boundary conditions. This is 
a subtle argument, because if we just start from the pe- 
riodic system, there seems to be no reason to expect the 
second term. In fact, this argument reinforces that the 
matrix element r x v cannot be used to represent the 
magnetization M in a setting with periodic boundary 
conditions, as the Mic term will be missing. 

(ii) The g — > calculation is not justified in the first 
place. As we know, even for free electrons, the wave 
function in an uniform magnetic field forms Landau lev- 
els, which are not pcrturbatively connected to the plane 
waves, in arbitrarily small magnetic fields. This is due to 
the fact that the perturbation is expanding in powers of 
A(q) in stead of B(q), where B(q) = q x A(q). If we do 
the perturbation formally anyway, in the limit g — > 0, but 
B remains finite, then the concern is that A(q) diverges. 
Indeed, in this setting we would find that the perturbed 
energy eigcnfunctions are not orthogonal to one another 
at any given order. Even though the correct formulas are 



recovered when the g — >■ limit is taken properly (prob- 
ably due to the fact that we are actually calculating the 
physical properties at B = 0), it is certainly desirable to 
have a more reliable derivation. 

(iii) The method of computing the response current 
to an adiabatic change of the Hamiltonian can be only 
applied to calculate the polarization. Magnetization 
change, for example, does not result in any bulk current 
fiow, and thus cannot be computed in any similar way. 
Another potential problem is that this method only cap- 
tures the change of polarization between the two systems. 
It is tempting, from the point of view of this method, 
to claim that only the difference of the polarization is 
physical. While one can always define the polarization 
of an atomic insulator to be zero and calculate between 
the interpolation of that and the state in question, it 
is not immediately obvious that any two state with dif- 
ferent polarization are measurably different when they 
are separately put with periodic boundary conditions. A 
derivation without referring to any other Hamiltonian is 
therefore desirable, as this directly shows that, for exam- 
ple, the polarization is an intrinsic property, independent 
of boundary condtions. 

With the issues mentioned above in mind, we would 
therefore like to develop a formalism, which explains 
and computes everything mentioned in a unified man- 
ner. In addition, since the integer quantum Hall effect 
and its higher-dimensional analogs are closely related to 
the quantities mentioned above and can be derived us- 
ing the Green's function techniques at finite momentum, 
we would like to propose a formalism utilizing Green's 
functions. In this paper, we provide such an unified 
formalism. In this formalism, we do not have to work 
with any boundary. We can also perturb in the uniform 
electromagnetic fields in a gauge-invariant way, without 
appealing to any finite momentum calculation. All the 
calculations are also done without changing the Hamil- 
tonian. The paper is organized as follows: in the sub- 
sequent section, we introduce the formalism and outline 
the procedure. In Sec. Ill, we show the derivation of 
various quantities in detail. 



II. FORMALISM 

What are polarization and magnetization? With 
boundaries, they can be defined as charges and cur- 
rents on the boundary; without boundaries, there has 
to be some inhomogeneity inside in order to observe the 
charges or currents. An alternative and more fundamen- 
tal definition is that the polarization (magnetization) is 
the coefficient for the term proportional to E(B) ,in the 
effective theory, when the electrons are integrated out. 
The boundary and the inhomogeneity charge or currents 
are then naturally derived when one solves the equa- 
tion of motion of the effective theory, which is just the 
Maxwell equations in our case. 

Therefore, our goal is to do the electronic part of the 
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path integral, in the presence of the uniform electric field 
and the magnetic field as a back ground, perturbatively in 
E and B. We then have P = -DF/dE, M = -dF/dB, 
and a) = ~d^F/dE'dB^, with F = fi-^ogZ the free 
energy. This at first seems rather straight forward, as a 
standard diagrammatic procedure is readily available to 
calculate perturbativc corrections to the partition func- 
tion. Our goal seems no more than a one-loop calcula- 
tion. It turns out not to be the case, however, when one 
looks carefully into the problem. The terms in the action 
of the effective theory we are after are total derivatives 
in terms of the electromagnetic gauge field. They are 
just zero in momentum space, where the standard pro- 
cedure is carried out. This also reflects the difficulties 
mentioned in the introduction, as either the operator r 
or r X V appear in the calculation of the polarization or 
the magnetization exactly due to the fact that E and B 
are spatial derivatives of the gauge potential. 

To overcome this problem, we have to calculate in po- 
sition space. Let us first deal with the magnetic field. 
While the wave function is not perturbative in the mag- 
netic field, as we will show below, the gauge-invariant 
part of the Green's function is. The Green's function 
for a single-particle Hamiltonian satisfies the following 
equation: 



H) xx' Qx'x" 



(1) 



H is the single-particle Hamiltonian and g is the Green's 
function of the electrons. Both are n by n matrices where 
n labels the orbitals and spins. The system couples to a 
small uniform magnetic field via the Peierls substitution: 



Fl^xx' — {H^ + H') 



A-dx 



(2) 



where Hq is the Hamiltonian in zero field, H' is some local 
perturbation that is proportional to B, e.g., the atomic 
diamagnetism, and A is the gauge potential. Since the 
correction to the Green's function as well as the free en- 
ergy from H' can be calculated in the standard way, we 
will set it to zero from now on. The line integral of the 
gauge potential follows a straight line from x to x' . In 
the following, we will use A^^' as the short-handed no- 
tation for A ■ dx. We also set e = h = 1 when there 
is no ambiguity. Using the idea in Ref.— , this equation 
can be solved perturbatively in B in the following way: 
we write 



9xx' — 9xx' ^ 



(3) 



and notice that it docs not change anything if we put 
g^^xx" along with the (5-function, we get 



"■ (4) 



Taking the exponential factor to the left-hand side, the 
three phases combine together, which gives the magnetic 



flux threading through the triangle formed by the three 
points x,x',x". Independent of the gauge, we therefore 
have 



Notice that this equation is now translationally invariant, 
and we can solve for g to flrst order in B by expanding 
the exponential and then Fourier transform, noting that 
X can be replaced by id/dk: 



9o 



1, iB-e-"- dg-^ dgo _ 
^ 2 dk'' dkb 



(6) 



where go = goik) = (a; - Ho{k)) ^] Ho{k) = 
Hofix cxp{~ikx) and g is in Fourier space. We there- 
fore get 



•^c^abc Qg-i 
5 -.90 + ^^.90^^ 



(7) 



Notice that g is gauge invariant. Once we have g{k), the 
Green's function is just the inverse Fourier transform of 
it times the phase factor e'^^="' . We therefore have the 
real-space Green's function in the presence of the uniform 
magnetic field. 

While the calculation is straightforward, to our knowl- 
edge Eq. ([7]) is a new result. In RefJ^, without sources 
other than the magnetic field which breaks time rever- 
sal symmetry, this first order term vanishes and all they 
have to do is to set g = go- In that case, all the effect of 
the magnetic field comes from the phase. 

We can extend the calculation to include the pertur- 
bative correction in the uniform electric field as well. 
We start from the defining equation which is the Fourier 
transform of Eq. ([l} : 

- H)(^x,x';t,t')9ix\ x"; t', t") = Sxx"Sit - t"). (8) 

Now we assume the coupling to the electric field comes 
from the space-time extension of the Peierls phase. Note 
that this procedure again does not include contributions 
from the response of the local orbitals to the electric field. 
We then use the same trick, define 



gxx'.w — gxx'.tt'G 



(9) 



where Axx',tt' is the line integral of the spacctime gauge 
field {~-V,A) on the straight line connecting the two 
points. Following a similar procedure, noticing that 



0% 



I, one can reach 



iE'' , dgo dgo ^ , .^.j^2n 
50-— (5o^-^.9o)+0(i? ). 



(10) 



This procedure can easily be carried to arbitarary order 
of both the electric field and the magnetic field. 

It is important to understand that Eq. ^ and Eq. © 
are just a way to factor out the gauge dependence of 
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the Green's function; it is not an approximation. The 
only approximation comes in when we Taylor-expand in 
powers of the flux threaded in the triangle formed by the 
three points. 

Let us be concrete and give a specific example. Sup- 
pose we have a tight-binding system with n orbitals sit- 
ting on each site. The i-th orbital is located at di from 
the lattice vector R. Now the Hamiltonian Hq is an n x n 
matrix in momentum space, and so is g. Notice that in 
deriving the formula, we have implicitly chosen the gauge 
such that -ffo(fc + G) — U'^Ho{k)U ; U is a. diagonal ma- 
trix with Uii = exp(iG • di)). The boundary condition 
is similar for g. Our formula is then a matrix equation 
for the n X n matrix g. It is important that our formula 
only works with this "twisted" boundary condition when 
there is a basis. 

We note that it has been shown earlier that the one- 
particle density matrix (OPDM) is also perturbative in 
the magnetic field, and can be calculated in a similar 
wayiii Many quantities we calculate below can also be 
calculated using the OPDM. One key difference is that 
the Green's function can also be perturbed in powers of 
the uniform electric field as we have shown above. Gom- 
bining the Berry's phase procedure as we will mention 
later on, the Green's functions formalism is thus a truly 
unified framework which can calculate perturbations of 
the uniform electromagnetic fields to arbitrary order, in- 
cluding the susceptibility and polarizability. The OPDM 
can always be derived from the Green's functions via 

At zero temperature, without the electric field, the free 
energy is just the expectation value of H. The path in- 
tegral can thus be performed by calculating the expecta- 
tion value of the Hamiltonian in a uniform magnetic field. 
Other perturbations in the presence of the field can be 
captured in the usual way, replacing fermion bilinears 
with the Green's functions. 

With a uniform electric field, the expectation value of 
the Hamiltonian is no longer the same as the free energy. 
We can understand this fact by taking the gauge V = 
(since our formulation is gauge independent.) In this 
gauge, the translational invariance in the time direction 
is lost, and one naturally does not expect any relation 
between the two quantities. One can directly see this 
by calculating the expectation value of the Hamiltonian 
in the presence of the electric field. We find that the 
expectation value of the Hamiltonian does not change 
with the electric field at first order. 

How do we calculate the path integral in the presence 
of the electric field then? The following observation pro- 
vides a hint: if we think about the imaginary-time path 
integral, the term P • E, unlike M • B, stays imaginary. 
Indeed as we discussed in Ref.~, the polarization 
is better thought of as a Berry's phase, instead of en- 
ergy, when the gauge winding number in the i-dircction 
is changed by one; i.e. 



Similarly, the extra Berry's phase in a magnetic field is 
related to the OMP by 



B crry 



(12) 



where <I>b is the total magnetic flux threading through 
the system. The 27r ambiguity of both quantities thus 
comes naturally. For the sake of completeness, let us 
reiterate the procedur o^^i^^ : 

Consider a system with periodic boundary conditions. 
To calculate accumulated the Berry's phase during a time 
when, for example, J A^dx increase by 27r, first we shall 
consider how the single particle wave function changes 
as we increase uniformly. The Bloch wave function is 
given by 



tpnkix) = Unk{x)e'' 



(13) 



where n is a band index; u„k{x) is periodic and satisfies 

((V - (fc + A,))' + F(x)) Unk{x) - EnkUnkix). (14) 

As we increase A^ uniformly to A^ + rj, the momentum 
k cannot change as it is fixed by the finite size L and 
the periodic boundary condition. On the other hand, 
following Eq. u„fc(x) changes as 



Unk{Ax + if) = Un(k-ri){Ax), 



(15) 



which is just a corresponding shift of k by —77. if 
77 = 27r/L, the system returns to its original state, but 
in a different gauge (i.e., with winding number different 
by one.) Notice that while u„kix) goes to the next ava- 
iable value on the left, the k in the exponential stays the 
same. The electronic wave function is therefore different 
from its starting state. Nevertheless, if we include the 
gauge field, the final state differs from the initial state 
by a large gauge transformation, and the Berry's phase 
accumulated in the process is well-defined. 

Now we are ready to calculate the accumulated Berry's 
phase of the band electrons under the process, where the 
winding of the gauge field in the x-direction is increased 
by one. Let us define A^ = J dxA^: 



"Berry 



f2TT{m + l) Q 

11 dA,(vI'e|-^|vI'e) 
/ 2-Km C'Ax 
^27r(m+l) ^ ^ Q 

oAx 



2izra 



ki .nGocc 



ki.n^occ ^ 



ki+2Tr/L 



ok^ 



/ dfc {Unk \ \Unk) 



(16) 



Berry 



-27rP,;. 



(11) 



\^ e) is the total electronic wave function; in the case we 
are interested it is just the Slater determinant of the oc- 
cupied electron wave functions at the wave vectors k^ al- 
lowed by the periodic boundary condition. In the second 
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equality, we wrote the derivative acting on the Slater de- 
terminent as a sum of derivatives acting on single particle 
wave functions. In the third equality we then plug in the 
dependence of the wave functions, and change variables 
to k. Whenever increases by 2tt, each Unk reaches the 
next allowed eigenstate to the left by the periodic bound- 
ary condition (without actually changing the momentum 
eigenvalue.) As we sum over all the integral of cigenstates 
at different allowed fc's, the whole Brillouin zone (BZ) is 
covered exactly once and we reach the fourth equality. 
We can read out the expression of the polarization using 
Eq. (HI]): 



px 



dk 

BZ 27r 



E 

nGocc 



Unkl 



_d_ 



\Unk) 



(17) 



This well-known result was derived in Ref.— i^^. 

How do wc express the Berry's phase in terms of the 
Green's functions? A naive thought would suggest that 
we cannot! Consider the "gauge transform" defined by 



(18) 



\Unk) exp(l0(fc)) \Unk) ■ 



Wc first observe that the Berry's phase, Eq. ([T6|) is 
not invariant under this transform and may change by 
integer multiples of 2tt. The Green's function, on the 
other hand, is clearly invariant under this gauge trans- 
form. It is therefore impossible to express the Berry's 
phase solely in terms of polynomials of the Green's func- 
tions. However, we then observe that the Berry's phase. 



Physical BZ 



Extended space 



Trivial system 




is hidden in the way we choose to extend the wave func- 
tions. The integrand on the surface can then be expressed 
in terms of the Green's functions, whose definition is also 
extended from the circle to the cylinder. The expression 
is as follows: 



Berry 



Uka) 



= i I d^ke'^9i (uka| dj |uka) 

JS 

= ]-[ d^ke^-'Tr' (ga,g-i9,g) ; (19) 
<J s 

dS is in the i-direction, which is the direction of the 
electric field, and S is the enclosed cylinder in the extra 
dimension. Bold quantities arc extended into the en- 
closed surface. We require that the other boundary we 
extend to docs not contribute to the Berry's phase. The 
trace with prime sums over all the bands (both occupied 
and unoccupied), and integrates over uj /2it as well as the 
other [d — 1) directions prependicular to E. Please see 
Appendix A for a derivation. 

This construction naturally separates the integral into 
two contributions, as discussed below: for the expression 
to be dependent only on the boundary variables, the in- 
tegrand on the surface has to be closed, in differential 
geometry terms. It is however, not always exact, in that 
integration over a closed surface does not always give you 
zero. A familiar example of an integrand being closed but 
not exact is VO on a circle, where 9 is the polar angle. It 
is locally a total derivative; nevertheless when you inte- 
grate it over the entire circle it gives you 27r. 

Wc can therefore separate the integral into two contri- 
butions. The first is exact and can be written as a total 
derivative in terms of Green's functions. It thus directly 
becomes a boundary integral via Stokes theorem. The 
remaining contribution cannot be written as a boundary 
integral in terms of the Green's functions. However, it 
has to be topological, meaning that it is invariant under 
smooth deformations which vanishes on the boundary. 
Topological integrands are of a specific form, as we will 
discuss in Appendix A. Specifically, in d spatial dimen- 
sions the topological term is in the form of the Wess- 
Zumino-Wittcn (WZW) actionii in d -I- 1 spacetime: 



FIG. 1. The black circle is the original physical space, with 
only the momentum direction along the electric field shown. 
The cylinder is the extended space, with the other end a trivial 
system. i,j are orthonormal basis on the extended space;i is 
along the direction of the eletric field, and j points along the 
extra dimension. 

expressed as a gauge-dependent loop integral in momen- 
tum space, can be cast as a gauge-invariant surface in- 
tegral via the Stoke's theorem. We therefore extend our 
system to one extra dimension in the momentum space, 
interpolating between the original system and a trivial 
system whose Berry's phase is taken to be zerojiii^ See 
Fig. 1 for an illustration. The gauge dependence, then. 



d+2 



■O'd + 2 



d 



dS 



UiUg^U-')^ (20) 



which is also defined in one extra dimensions. Our topo- 
logical term will be of the same mathematical form, with 
the spacetime x'^'^^ replaced by the physical momentum 
k'^ as well as the energy u, and the group element U{x) 
replaced by the Green's function go{k,uj). Further more, 
the coefficient in front of the topological term is deter- 
mined up to a sign by requiring that the value of the 
integral has a 2t: ambiguity with different extension to 
the extended space. This point was made in Rcfi^ for 
the case of time reversal invariance. 
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Let us summarize for the formalism: we have an ex- 
pression for the Green's function in the presence of an 
uniform magnetic field. Without the electric field, we 
can calculate the expectation value of the Hamiltonian to 
get the logarithm of the partition function. We can also 
calculate the charge and current responses. If we want 
to capture the terms linear in the electric field, however, 
we have to calculate a Berry's phase, which can only be 
expressed in terms of the Green's functions in one extra 
dimension. In the following section we shall see detailed 
calculations for all the quantities mentioned above. 

III. DERIVATIONS 

In this section, we show in detail how we apply the for- 
malism given in the previous section, to three different 
quantities: the charge response to the magnetic field, the 
magnetization, and the orbital magneto-polarizability 
(OMP). We also discuss the 6'-term, which is the isotropic 
part of the OMP, in higher dimensions. 

A. Charge response in a magnetic field. 



field, as our derivation implies. Without boundaries, the 
charge density can only change by adding or removing 
bulk states abruptly, even though the magnetic field is 
small and does not affect the energy gap in any apprecia- 
ble way. In a tight-binding model, the states then must 
"teleport" between the occupied and empty bands. On a 
torus the magnetic fiux has to be quantized; in the weak 
field limit, when we increase the magnetic flux by one fiux 
quantum, there will be exactly Ci states, "teleporting" 
from the unoccupied bands to the occupied bands, with 
out changing the energy gap in between. We have ver- 
ified this phenomenon with a numerical diagonalization 
of an insulating system with Ci 7^ 0. 

This peculiarity becomes more apparent when one 
compare the findings with the usual linear response 
derivation. There the quantized conductivity or den- 
sity change is derived from a bubble diagram at finite 
5, in the g — >■ limit. The density modulates in the 
same way as the magnetic field, which becomes uniform 
only at the limit. There are no such teleportations of 
the states between the bands; the electrons flows from 
patches with a positive magnetic field to patches with a 
negative magnetic field (if the Chcrn number is positive) 
and vice versa. 



In an integer quantum Hall system in two dimensions, 
the magnetic field is locked with the density. This is 
also true for Chern insulators, which has a non-vanishing 
transverse conductivity in zero field. The transverse con- 
ductivity is related to the field derivative of the density 
by the Streda formula a^y 
result as a sanity check. 

Starting from eqn. ([7]), it is straight forward to calcu- 
late the charge density in the magnetic field: 



^ . Here we first verify this 



Po 



(2^ 2^' tr yg, + —-g,^ — 
BCi 



2tt ' 



(21) 



Ci is the Thouless-Kohmoto-Nightingale-dcn Nijs 
(TKNN) index (Chern number) of the occupied bands. 
The capital trace in the first line implies summing over 
all the bands and integrating over all directions as well 
as energy with a facter of (27r)~^ each. The lowercased 
trace in the second line implies summing all the bands 
only. The phase factor in eq. ^ is absent due to the fact 
that the creation operator and the annihilation operator 
are at the same position. The last equality is derived in 
Appendix A. Eq. ([?!]) describes the locking of the density 
to the magnetic field. This "incompressibility" is funda- 
mental to quantum Hall physics and follows from charge 
conservation and Faraday's law when the magnetic field 
is adiabatically turned on. 

Nevertheless, it is somewhat intriguing to see this ef- 
fect survive even on a torus with an uniform magnetic 



B. Magnetization 

As sketched in the previous section, we calculate the 
orbital magnetization by computing the energy of the 
system in the presence of an uniform magnetic field, and 
the relation M = — d(iJ)/dB. We continue our deriva- 
tion in two dimensions, as the magnetization is inherently 
a two-dimensional phenomenon. We have 



E 

xx' 



En - 



Auj 



27r 



i tr {Hr^^ig^'x) 
d^k duj 



E 

xx' 



duj 



i tr (Hoxx'gx'x) 



(27r)2 2tt 



tr Hogo 



2tt 

dgo 
dk" dk'> 



(22) 



Plugging in go Y.m 1"" 



^■g— and noticing that 
the w-integral restricts the poles of the two go's to be on 
opposite sides, we find that the derivative on g^ can only 
act on the bra, and the derivative on g^^ can act on the 
bra or ket but not the energy, in order for the expression 
not to vanish. We can simplify the term linear in B to 
be 



(AH) 



Be 



ab 



d^k dw 



2 J (27r)2 2tt ^ 

Em{En — Era) 



[lO - E,n)iu} - En) 

iBe"^ r d^fc 



EmEn) (23) 
Um\daU„) {dbUn\Um) 



i2ny 



d{-E„iEn) 



X \E^\ ( \daUn) {dbU 
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In the derivation we have used relations such as 
{daUjn\un) = - {um\daUn) ■ Now we reexpress every- 
thing in terms of occupied bands only: 

^9{-EjnEn)\Em\ {Ura\daUn) {dbUn\Um) 
= ^ E„t {Um\daUn) {dbU 

n I '^m / 

nGocc,77iGcmp 

~En{u„\daUm) {dbUm\Un) 

= ^ {dbUn\H\daUn) - En' {Un'\daUn) {dbUn\Un>) 
n,n' Gocc 

-En{daUn\dbUn) + E^ {daUn\Un') {Un'\dbUn) 

= ^ {dbUn\H\daUn) - En {daUn\dbU„) . (24) 
nGocc 

Putting cq. back into the expression for the energy 
and dividing it by B, we recover the result derived in 
Ref^>^: 

M = ^ y ^ 1^ + EnldbUn) (25) 

nGocc 

Specifically, the first term is called Mlc and the second 
term is called Mjc in Refi^. One can also derive the same 
expression using the OPDMJ^. 

From our derivation, both terms come together from 
one simple expression in the bulk. If we have a boundary, 
the current flowing on the edge is the sum of the two, as 
is required by the Maxwell equations, J = V x M. One 
might be tempted to conclude that the two terms gen- 
erate separately measurable currents from the derivation 
in Refj^, but that is not the case. One can, however, 
decompose it into two contributions (different from but 
related to Mlc and Mjc) where the difference between 
the two terms can also be measured.— 

In addition, as also discussed in Ref.— , eq. (f25|) nat- 
urally shows that there are Ci gapless edge states for 
a Chern insulator with Chern number Ci in two spa- 
tial dimensions. Suppose we shift the chemical potential 
by A/i. Now the magnetization changes by CiA/i/27r, 
which implies the current on the edge changes by the 
same amount. If we think from the edge perspective, a 
chemical potential change of Afj, implies that the density 
of occupied edge states increases by CiAij,/2ttv, and car- 
ries additional current CiA/i/27r. The two observations 
would not have matched, had we only taken the contri- 
bution Mjc as the current flowing around the edge. 

On the other hand, we have the same orbital magne- 
tization formula for a system without boundary. This 
means that if the system is a Chern insulator, its mag- 
netization will depend on the chemical potential, even 
if there are no edges. This somewhat puzzling observa- 
tion actually comes from the density locking to the mag- 
netic field mentioned in the previous subsection. Since 
we define the magnetization as the derivative of the free 
energy with respect to B, with F ~ E ~ fiN, the free 
energy change with the magnetic field docs depend on fi. 



It is, however, unclear to us whether the "magnetization" 
defined this way is a measurable quantity on a periodic 
system. 

With the discussions above, wc therefore conclude that 
the magnetization of an insulator is a bulk quantity. The 
current flowing on the edge can always be derived via 
the equation of motion of the effective theory, i.e., the 
Maxwell equations. We also note that the orbital mag- 
netic susceptibility can be calculated by expanding the 
Green's function to second order in B. 



C. Electric Polarizability 

The calculation of the polarization is already covered 
in the previous section and in Reff^. While we can ex- 
press it in terms of the Green's function in one extra 
dimension as done in Eq. (|19l) . after integrating out lo 
and integrating back to the boundary, the result is just 
what wc start with. As a nontrivial example of using 
the formalism, here we calculate the polarizability, by 
considering the first order correction of the Green's func- 
tion in a uniform electric field. We start from Eq. (|19p . 
plug in Eq. (jlOp . and notice that there is an additional 
phase proportional to the electric field from contracting 
the three green's functions in real spacetime. (Notice that 
the time-dependent gauge potential does not break trans- 
lational symmetry in the spatial direction so Eq. (|19p still 
applies. The trace in the time direction, however, has to 
be carried out in real space; it is easier just to imagine ev- 
ery thing is done in real space and then converted back.) 
We thus have 

^<^Borry h + h] 

h = ^y d2ke'JTr'(^(goaago - da^o^o)d^^^^dj^o 

+ gOl9jgo ^l9j(gol9ago - 9agogo)^ 

h^-'-fJ^ d^ke'^^Tr' (^d^godad^g^'d.go^ . (26) 

Taking advantage of the relation di^g^^ — 1 and dujgo = 
— gQ, we can simplify the expression to 

A./'Lrry = ^ d^ke^^Tr' (^daid^godigodjgo) 

+ 2dt{godadjgo)^ 

Notice the only total derivative that does not vanish af- 
ter integration has to be along the extra dimension. As 
we see here, there is no topological contribution in the 
polarizability; the integral thus has to reduce to a to- 
tal derivative, and wc can just integrate it back to the 



8 



physical Brillouin zone. Using Eq. (fTTj) . wc get the po- 
larizabihty tensor 

which is the same expression as one would get using the 
usual perturbation theory at finite momentum q, and 
take the q ^ limit. 

D. Orbital Magneto-Polarizability 

Wc can calculate the OMP by considering the polariza- 
tion in a uniform magnetic field in three spatial dimen- 
sions. We start from Eq. ((T^ and plug in the Green's 
function in the presence of the magnetic field. Recall 
that i is the direction of the electric field, and j is the 
direction of the extra dimension: 

•/-Berry = ^d^ke'^Tr' (g^.g-^^.g) 

= 2n\'^L^TT'' {gd,g-'d,g) (29) 

Here g denotes the Green's function at a given momen- 
tum in the (i, surface. It is in general not translation- 
ally invariant in the remaining directions. Here we intro- 
duce the new notation Tr"^ for later convenience. Tr"^ is 
defined as integrating over momentum divided by (27r)^ 
in the direction within the boundary, integrating 

over aj/27r, and summing over positions, divided by L in 
the remaining directions if they are not translationally 
invariant; the lattice is replaced by integration over the 
momentum and divided by {2tt) otherwise. When E and 
B are not prependicular, we have to take a Landau gauge 
to make the Green's function translationally invariant in 
the direction of the electric field for the derivation in the 
previous section to work; it nevertheless does not affect 
the result, Eq. (P^]) . and our calculation below. 

Now we plug in Eq. ([7]). To first order in B, not only 
do we have g to first order, but we also have to con- 
sider that the three product of g's, contains three phases 
which sums to be the flux threading through the trian- 
gle. Indentical to what we did in Sec. H, we Fourier 
transform, Taylor-expanding the phase to first order in 
B. The result is 

^'/'Bcrry = ^1 + h] 

+ godiga^djigodago^dbgo)) 
h = -TT^iB-i^^-^-eyTr^ [dagMgo^d.ga) ; 

(30) 

Ii is from first-order terms in g and I2 is from the phase, 
go is the Green's function in zero magnetic field in the 
extended space. The trace now indicates integration in 
all the momentum directions as well as uj and divided by 



(27r) in each direction, as the translational invariance is 
restored. Ii can be rewritten as 

h = ^^iB-L\'^^-e'nr^{godag-^^dhg^d,g^^d,go 
+ {i^ a,j ^b)- ^■i{^^gQ^ago^^bgo)) (31) 

in which the last term in the second line can readily be 
integrated back to the physical momentum space. The 
first two terms are almost in the form of the topological 
terms we mentioned in Appendix A, but without com- 
plete antisymmetrization among the indices. 

Let us now look at l2- We would like to separate this 
term into a total-derivative and some remaining parts 
which give the topological term; to achieve this, we need 
to take advantage of the condition du:g^^ = 1. Note that 
this relation also implies d^jdago^ = and di^go = — go- 
By inserting di^g^^ at the end of the term and integrating 
the w-integral by parts, we find (from here on we omit 
the subscript of go to avoid cluttering in the equations): 

= Tr^(^g9,g-ig(afca,g-i)g9,g-ig 
+ g{dbdig~^)gdjg~^gdag~^g 
+ gdjg^^gdag^^g{dbd,g^^)g^ . (32) 

For convenience, we use the notation (abij) to stand for 
Tr'^{gdag~^gdbg~^gdig-^gdjg-^g), (a[bi]j) to stand 
for Tr'^(g9ag"^g(9&9ig~^)g9jg~^g), etc. The equa- 
tion above becomes 

(a[bi]j) + ([bi]ja) + (ja[bi]) =0. (33) 

Also, Eq. (|5T|) becomes 

(abij) + (ijab)-aj(iab); (34) 

We then integrate by parts twice on the second and the 
third term, noticing that a,b and i,j are separately an- 
tisymmetrized, to make them into the form of the first: 

([bi]ja) = (a[bi]j) + (ibja) - (bjia) + 9i(bja); (35) 
(ja[bi]) = (a[bi]j) + (jaib) - (jbai) - 9j(bai). (36) 

We therefore have 

- (a[bi]j) = J(-(bjia) - (jbai) + (ibja) + (jaib) 

+ ai(bja)-aj(bai)). (37) 
Now we can sum over all contributions, and get 

A<^j3,„y = 7:'^B''■Lh^'>^e^^ (^(abij) + (ijab) 

+ i((bija) + (jabi) + (ibja) + (jaib)) 

- 5j(iab) + i(ai(bja) - d^ihm))) (38) 
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Notice that the integrating-by-part triek in the lj- 
direction can also be apphed to expressions such as (abij) 
and (iab), and similarly we get 

(abij) + (bija) + (ijab) + (jabi) = (39) 
(aibj) + (ibja) = (40) 
(iab) + (abi) + (bia) = (41) 

Therefore, by writing 

(abij) + (ijab) = i ((abij) + (ijab)))-^ ((bija) + (jabi)) 

(42) 



A0Wy,wzw = ^7THB^Lh'^''^e^^ ( (abij) + (ijab) 



the "topological part" of (/)Bcrry is 

1 
3 



+(aijb) + (iabj) + (ibja) + (ajbi) j .(43) 

Notice that this term is totally antisymmetric in all the 
indices. This is expected, as it is topological only when 
all the indices are antisymmetrized. It is also of the form 
of the WZW action if we put d^g~^ into the expression. 

Since the direction of the magnetic field is prependic- 
ular to the {ah) plane, it has to be in the (ij) plane, for 
the topological term not to vanish. This implies that 
only the component of the magnetic field in the direction 
of the electric field contributes in the topological part. 
Gathering everything we finally have 



A0; 



Berry, WZW 



(44) 



The remaining part can be reorganized using Eq. ()4ip : 
2nH 



A(f>' 



Berry, 3d 



-$^e'^''=e*^(6»j(bai) + 6>j(iab)) (45) 



Integrating back to the physical momentum space, we 
have 



Berry, 3d 



-<^>'se'''"'{{iab) + {bai)). (46) 



Here {iab) stands for Tx{gdig~^gdag'^^gdbg^^g) with the 
trace summing over the energy as well as the physical 
momentum directions, with (27r)^^ in every direction. 
The difference of a factor of {2ti) in front comes from 
the different number of (27r)~^ in the definition of the 
traces Tr'^ and Tr, in four and three spatial dimensions 
respectively. Combining, we thus have our final answer 
in terms of the Green's functions: 



a.ij = (q!wzw + asd) 



i 



--^CafccdTr {gdag gdbg gdcg gddg g)%; 



a3dij ^ ~^<^ab]T^y{gdig ^gdag ^gdtg ^g-h.c). (47) 

Notice that in terms of the Green's functions, awzw can 
only be expressed with the extended dimension. Eq. (|47p 



generalizes the result in Refii to the generic time-reversal 
breaking cases, where the WZW integral can take con- 
tinuous values. Note that there is an additional term asd 
which is zero in the time reversal invariant case. 

To get the expression entirely in terms of variables in 
the physical momentum space, we have to expand the 
Green's functions explicitly in the eigenbasis then inte- 
grate it back to the physical momentum space. Taking 
advantage of the topological property of awzw and using 
Eq. (|63p in Appendix A, we can immediately know that 
the first term contains a part which can be expressed 
using the Berry's phase gauge field strength, and some 
other part which is a global total derivative, and invariant 
under the gauge transform defined by Eq. ([T8| . When in- 
tegrated back to the physical momentum space, the first 
part becomes the Chern-Simons term with the Berry's 
phase gauge field A^^^nn' = {unk\ - igfjr the re- 

maining part combined with a^^ gives the rest of the 
tensor aij as derived in Ref.— : 

o-tj = {acs + aQ)ij; 

1 /■ d'^fc 2 
acsi] = ^2^^] j j2TTf'^''^'''^^^'^''^^'^'' ^ * g AAA); 



•J \ / mP-cmTi Til 



+ C.C. , (48) 



mGemp,n£occ 

{d,n\m) {m\{daH,dbV}\n) 



where V — X^neocc I"") (""I projector to the occupied 

bands. The detail of the calculation is shown in Appendix 
B. 



E. ^-term in higher dimensions 

While it is straight forward to generalize the complete 
calculation in the previous section to obtain all compo- 
nents of the anolog of OMP in higher dimensions, the 
totally-antisymmetric part, i.e., the 0-term, is especially 
easy to compute. Here as an illustration, we calculate the 
coefficient 6^^ of the 0-term in {5+l)D as defined below: 



^5d 

3847r3 ' 



(49) 



First we find the Green's function to second order of the 
magnetic field. In higher dimensions, the term in the ex- 



ponential in Eq. ([5]) becomes [iFab{x' — x)°-{x'' 



72]. 



Since we are only interested in contributions with all the 
indices antisymmetrized, only first derivatives will con- 
tribute. When we Taylor-expand Eq. ([5]) to second order 
and Fourier transform, we get 



.92 = —^{FabFcd + FadFbc + FacFdb) 

90 {dago^)go{dbgo^)9o {dcgo ^ )go {ddgo ^)9o 



(50) 



where 52 is the second order term of g, and the indices 
run through all five spatial directions. The (...) vanishes 
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when wc antisymmetrize all the indices. Plugging into 
Eq. ([TO]) , we then have 



3^2 
2 

(ijabcd) 



^Bcrry — ^e.'^'' (Fab Fed + FadFbc + FacFdb) 



(51) 



where we have used the abbreviated notation introduced 
in the previous section. We still only need to keep track 
of the parts which do not vanish after antisymmctrizing 
all the indices. Antisymmetrizing, noticing that in each 
direction (say E^, By^, B^v) summing over indices gives 
a factor of 8, we then have 



Berry 



abode f , 



60 



abcdef)$^$|, 



(52) 



and $^ arc the two magnetic fluxes threading 
through the four directions prcpcndicular to the electric 
field. Notice that from the definition Eq. p9|) . 9r.ri is ex- 
actly the Berry's phase when the flux threading through 
each direction equals 27r, we then have 



75d 



15 



1 



19207r 



(53) 



Ip is defined in Eq. (|60| . This is the higher-dimensional 
analog of the trace of awzw, which includes both the 
second Chern-Simons term A /\ F /\ F, and some other 
inter-gap contributions. To rccxpress entirely in 
terms of Bloch wave functions and energies in the phys- 
ical Brillouin zone, however, is rather tedious, and we 
shall not do it here. 



IV. DISCUSSION 

In this paper, we provide a formalism to integrate 
out the electrons with external uniform electromagnetic 
fields. This formalism provides a unified and system- 
atic way to calculate quantum-Hall type responses, po- 
larization, polarizability, orbital magnetization, orbital 
magnetic susceptibility, and OMP. From the perspective 
of the formalism, all of the quantities mentioned are of 
bulk nature, and all calculations can be done with pe- 
riodic boundary conditions. The existence of the edge 
current or response can be derived from the equation of 
motion of the resulting effective theory, which is defined 
in the bulk, independent of the boundary conditions. 

In our formulation, one key insight is that the linear 
term in the magnetic field is an energy, whereas the lin- 
ear term in the electric field is a Berry's phase. This 
explains why the polarization and the OMP are defined 
only modulo 2tt in certain units, whereas the magnetiza- 
tion is always rigorously defined. 

The wave function under an uniform magnetic field is 
nonperturbative; nevertheless, the gauge invariant part 
of the Green's function is perturbative, even strictly at 



(7 = 0. In Eq. ([7]), the only expansion parameter is 
the flux enclosed in the triangle; inside an insulator with 
finite range correlation functions, the expansion is con- 
trolled. 

When the 6'-term is first discovered in the 3D topo- 
logical insulators, it was shown for time-reversal invari- 
ant systems, 6 is given by the Chern-Simons term of the 
Berry's phase gauge field defined in momentum space. 
It is not until much later that extra contributions which 
depend on the inter-gap matrix elements as well as the 
gap size are discovered for the general case without time 
reversal invariance, along with the off-diagonal compo- 
nents. From our calculation, two terms come together 
naturally, and are most simply expressed as an integral 
of the WZW term with Green's functions in the extended 
space. The same conclusion holds for the higher dimen- 
sional 9 term as well. 

However, we have to note that it is not easy to con- 
vert the expression in the extended space back to the 
Bloch wave functions and energies in the physical Bril- 
louin zone. For this purpose, the density matrix per- 
turbation theory formalisnii^ seems to be more useful. 
Nevertheless, from those methods it is harder to obtain 
the Chern-Simons term; it is also not as straight-forward 
to generalize to higher dimensions. Our formula Eq. (|47p 
and the higher dimensional generalization Eq. (|53p thus 
complement the other methods and offer a better con- 
ceptual understanding. 

There is an important difference between the deriva- 
tion of the 0-tcrm using our approach and the dimen- 
sional reduction procedure used in Refs.— i^. Even though 
they give the same result with time reversal symmetry, 
the latter is not readily generalizable to the general case 
without time reversal symmetry. In the dimensional re- 
duction procedure, one thinks of the system in one higher 
dimension as a collection of systems in the physical di- 
mensions with different 9. The sum of the 0-term of the 
systems in the physical dimensions gives rise to the quan- 
tized transverse response in one higher dimension: 



0'2d 



dk 89 



2d-l 



2tt dk 



(54) 



a2d is the coefficient of the quantized transverse response 
in 2d spatial dimensions, and k is along the direction of 
the extra dimension. For a2d to be nonzero, 02d-i has 
to be gauge dependent with the gauge transform defined 
in Eq. (fT8|). From the expression of a2d, all one can de- 
duce is that 92d-i must contain the {2d — l)-dimensional 
analog of cither a^zw or acs , as well as some arbitrary 
quantity which is gauge invariant. In three spatial di- 
mensions, 9 is odd under time reversal symmetry; all the 
gauge invariant terms then are required to be zero when 
the symmetry is present. This is the only situation in 
which ttwzw and acs, given in Eq. (|47|) and Eq. (|48|) 
respectively, are equal, and the derivation in Ref.— un- 
ambiguously determines the answer. 
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APPENDIX A: TOPOLOGICAL COMBINATION 
OF GREEN'S FUNCTIONS 



tice here we did not assume any particular form of that 
is, / is topological under arbitrary smooth deformation 
of g. 

We can construct similar topologically-invariant inte- 
grals in higher dimensions. In fact. 



2d+l 



In this section, we first discuss what kind of combina- 
tions of Green's functions are topological, then we pro- 
vide some formulas we have used in the main text. 

When we say an integral is topological, we mean that 
the value of the integral is independent of any smooth 

deformations that leaves the integrand on the boundary , , r , ^ -^^ i 

J 1 J.1 r n • • J. 1 • i we have an even number oi qc/n o multiplied together 

unaltered. I'or example, the following integral is topo- . . ' ^ , ^ , 

with their indices antisymmetrized, they become a total 

derivative: 



is topological. To prove this, we notice first that when 



(55) 



61 



logical: 

I = J duj tT{gd^g^^) 

duj tT{Sgd^g~'^ + gd^jSg'^) 
duj iT{6gdu:g~^ + d^{g5g^^) - d^gSg^^) 
dcj tr{Sgd^g-'^ + d^{gSg^'^) + d^gg^^Sgg'^) 
dcotiidUgSg-')); (56) 



2d 



d 

e^^'^^-'^^-^{~irY[{da,^_,g){da,^g-') 



(58) 



J=2 



we have used 6{gg ^) = duj{gg "'^) = in the last two Now we consider a general deformation: 
I 

» 2d+l 

81^" = (2d+ 1) / d^d^-^fc e'^i'^-'^-+Hr((5(gaai.g~i) [] i9da,g~^)) 

i=2 

/2d+l 
diod^^k e''^<'-----+Hv{{5gda,g-^ +daA95g-') + da,gg-^5gg-') \{ igda,g-')) 
i=2 

/2d+l 2d+l 
du^d'^'k e^^^--^^^+Hii{6gda,g-' ~ gda,g-'6gg-^) J] {gda^g-') + daAgSg-' H (dda.d-'))) 
1=2 1=2 

» 2d+l 

= (2d+ 1) / dc^d^'^fc e'^i'^— -+Hr(Oa,(5,5g-i J] igda,g-'))); (59) 



the second-to-last equality follows from the fact that the 
product term is already a total derivative and with an- 
tisymmetricity of the indices, the product of two to- 
tal derivatives result in another total derivative. The 
last equality comes from the cyclic property of the trace 
and the first two terms cancel each other. We therefore 
have showed that l2d is topological for any d. Notice 
that even when g satisfies a twisted boundary condition 
g{k + G) = Wg{k)U on a closed manifold, the variation 
still vanishes. 



In the discussion we have in the main text, how- 
ever, we have concentrated on non-interacting systems, 
where the Green's function is taken to be in the form of 
(u! — Hk ± We thus are more interested in integrals 

of Green's functions that arc topological subject only to 
any smooth deformation of Hk, instead of an arbitrary 
deformation of g. Fortunately, since d^g~^ ~ 1 with the 
non-interacting Green's function, we can directly trans- 
late the topological combination above to combinations 
that is invariant under only the deformation of ■ The 
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invariant reads 



T2d 



2d 



M(R(9da,g-'))g). (60) 



i=i 



The indices now only run through all the spatial direc- 
tions. 

Next we show the derivation of Eq. (fT9)) . Since the 
Berry's phase is invariant under small deformations of 
the Green's function extended into the surface, the inte- 
gral has to be topological. By dimensional counting we 
immediately see that Ip is a possible candidate. The rest 
of the task is to find the constant in front, as well as any 
possible total derivatives. 



d^fce*^ 



d^fc 



cVk 



Tr {gd,g~'djg] 
'^'^ "tr {gdig' 



2tt 

— ( 

27r 



E 

mn 

E 

mn 



'd,g) 

\diHk\Un) {djUn\Urn) 



. — Em.) 



\diUn) {djUr 



\Ur,i) 



d^ke 



E 



En) 



m£cmp,?i£occ 



E 

nGocc 



d'ke'^ V i{d,Un\djUn) 



\diUn) 



(61) 



As we can see, there are no extra total derivative terms. 
If we integrate this term on a torus, as we did in Eq. ([2T|). 
we get 



Similarly, we can relate Ip to C2, the second Chern 
number if integrated on a four dimensional manifold 
without boundaries. In Refj^, they use dispersionless 
bands to derive the ratio betwenn if, and C2. Here we 
just quote their result (with an extra factor of (— i) since 
they are using the imaginary time Green's function): 



C2 



1 



327r2 

i 



d^ke 



ijki 



kl] 



r4 

487r2^^' 



with 



-|- z fli, a,- 



(63) 



(64) 



is the Berry's phase gauge field strength. Rather remark- 
ably, the two integral no longer agree with each other 
when integrated on a manifold with a boundary. In this 
case, they differ by a globally defined total derivative. 
APPENDIX B: FROM THE GREEN'S 
FUNCTION EXPRESSION OF THE OMP TO 
THE EXPRESSION IN ENERGY EIGENBASIS 



Before we start to evaluate equation (|48p . we empha- 
size again that our distinction between atopo and aoxact 
is different from the distinction between acs and uq in 
Refjii. Specifically, aexact is traceless, as can be seen 
from Eq. (|4ip . whereas ac in general is not. With this 
in mind, let us start from aoxact ^ 



--eabj{{iab) 



(65) 



BZ 



d^fce^^Tr {gd^g-^djg) = / d^fce 



BZ 

27rCi. 



As done in the main text, we use {iah) in short for 
Tr{gdig^^gdag^^gdbg^^g)- The capital trace again de- 
i {diUn\djUn) notes tracing over all the bands and integrating over all 
n6occ the momentum and energy divided by (27r) for each di- 

(62)rection. Expanding in the energy eigenbasis, we get 



, f d^k duj ^ {m\d,H\n) (n\daH\i) (i\dbH\m) 



(2^)3 27r ^ {lU - Ern)H^ - En)iLO - Ei) 



d^/fc / 1 

sgn{En)e{-EnEe)0{-E„E„,)- (9»m|n) {dan\e) {£\dbH\m) 



(27r)3 V En - E,n 



sgn{Ei)ei-EeE„,)ei-EiEn)—^——{m\d,H\7i) {dan\£) {dbi\m) 
sgn{E,n)e{-EnEm)e{-E,Em) ^ (a,m|n) {n\daH\t) {dd\m) 
sgn{Em)O{-EnEm)0i-E,Em)—^-—{dMn) {n\daH\i) (a,£|m) 
* / E ( e -E {dan\m') (m' \dbH\m) 

rii,7n' Gcmp;?i,n'Eocc 
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{m\diH\m) {dam'\n) {dbn'\m) -— — {dim\n) {n\daH\n) {dbn'\m) 

— — — {dim\n) {n\daH\n) {dbn\m) ) — (occ ^ emp); (66) 

we have used \n) to denote \unk) to avoid cluttering the expression. 

Since ctwzw is isotropic, the ofF-diagonal tensor element comes entirely from a^d- Without loss of generality, let us 
look at a„ 



lyz- 

i 



ctyz = --e°''^((?/a6) -c.c.~ 



m,m' Gcmp;7i,n' gocc 



= --^[[{.yxy) - {yyx)) - c.c.y (67) 

In the last equality of Eq. (|66)) . the first and the last term arc already in the form as derived in Refj^. Let us look 
at the two terms in the middle as we first plug in i ~ y, a = x, b ~ y: 

(— — ^——{m\dyH\m') {dxm'\n') {dyn'\m) — (occ -H- emp) ) 

m,m Ecmp;n,n tEocc 

(— — ^ {dyn'\m) {m\dyH\m') {dxm'\n') 

mp;n,n tEocc 

^——{dym\n) {n\dyH\n') (9,n'|m) (68) 

(— (9ym|n) {n\dxH\n') {dyn'\m) - (occ o cmp) ) 

( " ^^y^'\'^'> (9ym\n) {n \dxH\ n') 

.n'Gocc 

{dym\n') {dyn'\m') {m' \dxH\m) ). (69) 



m,m' Gcmp;n,n'Gocc 



ni,m,' Ecmp;n,n'Gocc 

1 



En' — Em 



Similarly, when we plug in i = y, a = y, b = x and sum with its complex conjugate, we find that every term cancels 
out: 



\ F j_ F {m\dyH\m') {dym'\n') {dxn'\m) - ^ {dym\n) {n\dyH\n') {dxn'\m) 

(occ o- emp) ) + c.c. = 0. 



m,m'^cmp;n,n'^occ 



(70) 

Summing over, we then have 



( {yxy) - {yyx)) - c.c. 

— — -r — — {dym\n) {dan\m') {m' \dbH\m) - (occ ^ emp) - c.c; (71) 

m,m (Ecmp;n,n Eocc 

here a and b run through only x and y. We finally have 

"''^"2/ (2^ ^ _^ (a^m|n) {dan\m') {m'\dbH\m) - (occ o cmp)j + c.c. (72) 

which is identical to the expression in Refj^i. 

I 



As for the diagonal components, we have to expand the isotropic term awzw The calculation is similar to 
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Eq. (j66p . except that we have to rearrange terms into 
total derivatives and then integrate back to the physical 
momentum space. We currently are not aware of any 
special trick to automatically rearrange the terms into 
total derivatives other than the J- A J- part, which cor- 
responds to acs in Ref.— ; nevertheless, we can use the 
Stokes theorem to convert the difference between aQ in 
Refiii and to the extended space and verify that it 



agrees with the remaining parts in a^zw- 

1 f d^k , , {dcn\m) 

— (occ <-> emp) 

+C.C. — (occ o emp). (73) 
We now start to expand awzw^ 



bed ^ 5,„^ „-i ^ „-i ^ „-i .0 „-l„^ _ /• d^fc /• dc. ^ {^ \daH\ j) {j \dbH\ k) {k \d,H\ i) (£ \ddH\ i) 



d^fc 



{J(ik,j£) + (-^(ij./c^) - c.c.) + J(ijfc£) + + {J{j,ike) - c.c.))) - (occ o emp). 



(27r) 

(74) 

We separate the contributions by different pole placements. {ik,j£) for example implies that Ei and Ej are of the 
same sign and are opposite to Ei; and Ef. Here we illustrate the calculation by explicitly computing the contribution 
J{ijM)'^ the calculation for the other contributions are similar and we shall just list the answer. To avoid cluttering 
the expression we use i to denote Ei as well when there is no ambiguity. 

-^.fe,,.)-^|_,^^gn(j)(^ (,_,)2(^-_fc)(^-_^) + (^_,)2(^_j-)(^_fc) y ^^a^\J) {d,j\k) {dM) {ddV) 

= E -Bgn(j)( [':')^|':;j + ^z^z'^^ Y'' ^^m^) (^^^^n) 



{ij.ki) 



^ sgn(j) _ + (j o £)]e"^^'^ (9fcj|fc) {d^) {ddV) 

J2 sgn(j)[^-i]e"^^'^(5„i|.7) {d,3\k) (8^) (dam + c.c. 
ij,ki) \3 ^ / 

J2 sgn(jV ^~^ + p' -^V""-'^ {dat\j) {dbj\k) {d,k\£) (ddli) +C.C. 
J2 sgn(j)f^ + ^-f)e'^'"=''(aa*U) {dbj\k) {dck\£) (3,^^) + c.c. 

E 3\ ,,,, ^^^^^^^ ^^^^^1^,^ {d,n'\m') {ddm'\n) + c.c. 

Tn,m'Gcmp,n,7i'Gocc 

(occ o emp). (75) 



To get to the first line, we have used 

{l \daH\ j) = daEAj + [E, - E,) (dallj) . (76) 

In J[ik,ji) the first term on the right-hand-side does not 
contribute. In other parts however that term is impor- 
tant. In the fourth equality, we have taken advantage 
of the fact that exchanging j and £ gives the complex 



conjugate of the whole expression. We symmetrize the 
second fraction in the third equality with its complex con- 
jugate to get i. In the second-to- last equality, the two 
fractions become each other with a negative sign when 
we exchange the occupied and empty states, therefore 
reaching the last equality. 

The other contributions are given in the following: 
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^{Ji^J,M) - = E 2 f ^"' j?"' + e""'' (^'^'^l"') (^''"'l^^) (9cm|m'} (9,,7i'|n) + c.c. 

n,n Gocc,7n,m^cmp 

+ E 2 f -^-)^-;^» ) .'^^c^ (5,^1^^) (a,n|m) +c.c. 

nfi^occ;m,m^cmp 

— (occ -O- emp); (77) 



-*J(fc.,,£)= E f ^^" (^^"V> (9cm|m') {d,m'\n) 

71,71' ^n" Gocc;m^cmp 

+ E f ^^7r ) 6°"-'^ (9fcm|mO (a.m» (9,n|m) + c.c. 



nGocc;m,m'gcmp 

(occ o cmp); (78) 



-^%.fc^)= E _( (-^» (Jr-Il)2^"''0 ^^^"''""^ ^^^""'""'^ ^^""^''"^ 

n,n',n" Gocc;m^cmp 

E f ^;"^"" - ^""^"" V""'" (^^"1"') (^'^"'V) + ^-^^ 



n,n',n" Gocc,mGGmp 



E 2 f ^^7;^ A \ ^abc, ^g^^i^,^ ^g^^i^^ ^^^^^ 

^ \ l-^m ^71 ) / 

nGocc;m,m' ^Ginp 



nGocc;m,m' Gcmp 

— (occ emp); (79) 
- C.C.) = E 1^ - E'-E„ \^,^, ^^^^1^^ ^^^^1^,^ ^ 



E ( ^"(g-' ^^) \a,cd ^ ^^^^ 



nGocc;m,m'Gcinp 



n£occ;7n,m' t^cmp 

(occ o cmp); (80) 



r 



notice that every term (that is not canceled) come with 
its complex conjugate. The term in J^ik.jt) without en- 
ergy dependence is the TAJ- term since 

^ab' = E (^-"1"^) idbm\n'} - (a O b). (81) 

mf^ernp 

The remaining terms sum up to be the following total 



derivative, as one can verify by taking the derivative on 
every term in the bracket in the following: 



{ddn\m) 



En — Em , 

(82) 

By comparing this expression to Eq. (|73p. we finally re- 
cover Eq. 
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